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Summary 


lull  solution  of  plane  wave  diffraction  by  a  surfa¬ 
ce  impedance  discontinuity  on  a  plane  is  presented.  The 
solution  is  given  for  both  normal  and  oblique  incidence, 
the  diffraction  integral  being  evaluated  asymptotically 
in  a  closed  f&rm.  This  solution  extends  previous  results 
of  Maliuzhinets ,  valid  only  for  normal  incidence,  and 
singular  at  reflection  boundaries. 

For  practical  application  to  radiowave  propagation  over 
the  Earth,  the  field  close  to  the  surface  is  of  interest. 
This  is  represented  as  incident  .reflected  .surface  and 
extra  term  field,  this  last  characteristic  of  the  surface 
impedance  discontinuity.  In  addition,  a  cylindrical  wave, 
scattered  by  the  discontinuity  is  present. 

The  f  ield  produced  by  prescribed  sources  can  be  synthes¬ 
ized  using  a  plane  wave  expansion  and  then  applying  results 
of  the  scattering  analysis. 


key  words  :  radiowave  propagation,  scattering 
asymptotic  evaluation,  plane 
wave  expansion. 
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1 .  Introduction 

Groundwave  propagation  is  one  of  the  oldest  studied  topics 
in  applied  electromagnetics . Several  models  of  the  propaga¬ 
tion  medium  have  been  considered,  e.g.,  flat  or  spherical 
Earth , homogeneous  ,  stratified  or  anisotropic  ground,  corrug¬ 
ated  surfaces , single  or  mixed  path  propagation.  An  excellent 
summary  of  problems  and  modern  solutions  to  this  subject  is 
given  in  [ 1 J ,  which  contains  also  a  large  number  of  references 

A  canonical  problem  in  mixed-path  propagation  is  that  of 
an  abrupt  discontinuity  in  ground  parameters  along  a  straight 
line,  as  depicted  in  Fig.l  :  the  ground  exhibits  different 
properties,  e.g.,  conductivity  and/or  permittivity,  for 
>  >  0  and  y  <  0.  The  canonical  problem  is  the  follow¬ 
ing  :  for  a  given  incident  plane  wave,  compute  the  field  at 
point  P.  If  the  plane  wave  solution  is  needed  to  synthesize 
the  field  produced  by  localized  sources,  the  angle  <J>  should 
be  allowed  to  be  complex,  and  the  incidence  should  not 
be  restricted  to  the  normal  case  (two-dimensional  problem) . 

A  central  role  in  the  solution  to  this  problem  is-played 
by  the  surface  impedance  concept,  i.e.,  by  definition,  the 
ratio  af  the  tangential  components  Et  and  of  the  fields 
at  the  i  interface  x  ■  0.  For  the  problem  of  an  incident 
plane  wave  with  a  fixed  angle  <j>  ,  the  use  of  an  impedance 
type  boundary  condition  is  an  exSct  approach.  However  ,  when 
using  plane  wave  expansions  (as  in  the  solution  of  the  just 
mentioned  canonical  problem)  the  use  of  an  angle  in4epend- 
ent  surface  impedance  implies  an  approximation  to  the  exact 
boundary  conditions. 

A  discussion  about  the  validity  of  the  surface  impedance 
concept  is  given  in  [2^  and,  more  recehtly,  in  (jf] .  We  want 
also  to  mention  a  few  'of :  the  excellent  agreements  between 
experimental  and  theoretical  results  using  impedance  type 
boundary  conditions:  see  [4,53  for  the  case  of  scattering 
by  a  half-plane  with  two  face  impedances,  and  [6-8]  .for  the 
case  of  laboratory  models  of  groundwave  propagation,  Add 
itional  results  are  listed  in  [l]] .  Accordingly,  the  ..use  of 
impedaiy:e  type  boundary  conditions  seems  quite  adequate  to 
the  case  at  hand. 

The  solution  to  the  canonical  problem  depicted  in  ,Fig.l 
is  kn^wn  essentially  for  the  case  of  normal  incidence 

’  I 

'  j 

% 

f 
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[9-17]  ;  when  oblique  incidence  is  considered,  only  an 

approximate  matching  of  tangential  components  of  fields 
at  the  impedance  discontinuity  line  is  provided.  We  note 
that  the  problem  of  scattering,  under  normal  incidence, by 
a  wedge  of  arbitrary  angle  with  two  different  face  imped¬ 
ances  has  been  rigorously  solved  since  1958  by  Maliuzhin- 
ets  [2(3]  .  When  the  wedge  angle  is  made  equal  to  n  ,  the 
problem  of  Fig.l  is  recovered.  Maliuzhinets  solution  has 
been  recently  extended  to  ,  the  case  of  oblique  incidence, 
for  the  half-plane  case  [2lJ  ;  we  are  using  essentially 
tiie  same  technique  to  construct  the  solution  to  the  problem 
of  Fig.l  for  the  case  of  oblique  incidence. 

Before  considering  the  extension  of  Maliuzhinets'  solu¬ 
tion  and  its  modifications,  appropriate  to  groundwave  pro¬ 
pagation  problems,  it  is  worth  remembering  that  an  altern¬ 
ative  convenient  way  to  study  mixed-path  propagation  pro¬ 
blems  is  obtained  by  using  the  compensation  theorem  [22]]. 

This  procedure  has  been  originally  exploited  by  Wait  [23- 
24]  and  then  applied  to  oblique  propagation  of  groundwaves 
across  a  coastline  [25-27]].  Extensive  theoretical  and  expe¬ 
rimental  work,  using  essentially  this  procedure  is  referred 
to  in  Q  ].  Use  of  the  compensation  theorem  is  very  attractive 
since  it  allows  you  to  consider  not  only  the  simple  canonical 
problem  of  Fig.l,  but,  in  principle,  any  geometry  of  sur¬ 
face  discontinuity,  therefore  rendering  the  procedure  very 
powerful  from  the  application  viewpoint.  However,  the  result¬ 
ing  integral  equation  should  be  solved  numerically,  even 
when  some  simplifying  approximations  are  made.  Accordingly 
the  rigorous  solution  of  the  canonical  problem, e.g.  >  the  one 
depicted  in  Fig.l, is  worth  a  physical  description 

of  the  refraction  phenomena  at  the  impedance  discontinuity. 

in  the  following  we  will  present  this  solution  for  both 
normal  and  oblique  incidence,  using  an  uniform  asymptotic 
evaluation  of  the  diffraction  integral  in  order  to  obtain 
a  field  description  bounded  everywhere.  This  is  not  a  tri¬ 
vial  problem  when  the  observation  point  is  close  to  the 
surface,  and  the  incident  wave  is  at  grazing  angle,  since 
as  many  as  three  poles  cluster  around  the  "saddle  point". 

The  difficulty  has  been  overcome  using  an  "ad  hoc"  modified 
version  of  Bleistein  's  procedure  [3l]  ,  thus  obtaining 
pble  type  contributions  in  addition  to  the  usual  saddle- 
point  contributions  (cylindrical  wave  emanating  frofo  the 
surface  impedance  discontinuity) .  The  pole-type  contribution 
consists  of  the  incident,  reflected,  surface  and  an  extra- 
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term  field,  this  last  characteristic  of  the  surface  imped¬ 
ance  discontinuity. 

The  field  produced  by  prescribed  sources  can  be  synthesyzed 
using  a  plane  wave  expansion  and  then  applying  results  of 
the  scattering  analysis.  Evaluation  of  the  resulting  integral 
is  made  difficult  by  the  existence  of  a  rapidly  varying  Fre¬ 
snel  integral  in  the  integrand 


2 .  Scattering  under  normal  incidence.  Application  of 
Maliuzhi nets' theory 

Let  us  consider  the  geometry  depicted  in  Fig.l,  wherein  a 
plane  wave  is  normally  incident  at  a  (possibly  complex) 
angle  <p  on  a  two  region  half-space.  No  z-variation  of 
the  fielS  is  assumed  (two-dimensional  problem).  Accordingly, 
we  can  consider  two  (decoupled)  types  of  field  polari zation- 
the  E-type,  wherein  the  only  component  of  the  electric  field 
is  E  (p,4>),  and  the  H-type,  wherein  the  only  component  of 
the  magnetic  field  is  H  (p,4>).  In  order  to  siinplyfy  the  no¬ 
tation  we  will  use  the  symbol  F  for  the  field  assuming  Fe= 
=Ez  ,  and  F^ *  and  using  the  suscripts  only  whenever 
necessary. 

We  will  assume  the  boundary  conditions  at  <J>  =  *tt/2  to  be 
of  impedance  type,  e.g., 

1  3F  - 

-  liksm0  F  =  0  at  ♦  =  *  tt/2  (2.1) 

where : 

±  .  ±  ± 

sin  0g  -  C./Z-  ,  sin  6^  =  t/Z  (2.2) 

k  *  w/ey  and  is  the  free-space  propagation  constant, 

C  =  /y/e  is  the  free  space  intrinsic  impedance; 

Z*  are  the  surface  impedances  at  $s±tt/2  respectively, 
and  a  time  dependance  exp  (-iwt)  is  understood  and  sup¬ 
pressed. 

Let  us  now  express  the  field  as  a  superposition  of  plane 
waves,  hence  :  0 

F^P,<^  =  2  tt  i  \  F(p+4>)exp(-ikpcosa)  da  (2.3) 

r 

wherein  F  is  the  two-branch  Sommerfeld  contour  of  inte¬ 
gration,  as  depicted  ifi  JUg  2.  Forcing  boundary  conditions 
(2.1)  in  (2.3)  gives,  upon  integration  by  parts, 
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j  ( sirm  t  In  6  )F(a  *  tt/2)  exp  (-ikpcosa)  du  =0  (2.4) 

'Y 

A  necessary  and  sufficient  condition  for  the  integral 
(2.4)  to  vanish  is  the  integrand  to  be  an  even  function  of 
a,  hence  : 

±  ± 

(sina*  sine  )F(  a  1  tt/2)  =  (-sina-sin  0  )F(-a  -  tt/ 2 )  (2.5) 


due  to  the  symmetry  of  the  integration  contour  r  .  tq.  (2.S 
can  be  solved  by  means  of  the  method  of  logarithmic  deriva¬ 
tive  and  Fourier  transform.  Only  the  general  solution  is 
quoted  hereafter,  under  the  assumption  that  a  single  pole 

at  a=0  does  exist  in  the  strip  -  tt / 2 <  a'  <  n/2 
o 


F(a)  = 


sir.  a -sin  0  n=0 

o 


qp 

+  X  Bnsinna 


0(a) 


(  2  .  ^  ' 


wherein  A  and  B  are  arbitrary  constants  and  the  function 
n 

0(a)  is  introduced  and  discussed  in  Appendix  A. 

As  discussed  in  [29}  ,  the  value  of  the  field  F  for  p  — 0 
is  related  to  the  asymptotic  behaviour  of  (2.6)  as|a"|  — ►  cd 
In  particular  ,  F(p,<|>)  is  regular  for  p  — >  0  if  F(a)  is 
bounded  for  (a"  |  - »  cd  .  Since  from  results  of  Appendix  A 

.  0(  a  )  %  C  (  exp  ( |  a”  |  )  )  (2.7) 


as 


CD 


F(A>  =  F 


it  follows  that  =  0,  thus  [20J 

COS  0  ,  ,  . 

o  0(a) 


(2.8) 


sina  -  sin  0Q  0(0Q) 

wherein  F1  is  a  constant  (corresponding  to  the  intensity 
of  the  incident  field  at  px0  ,  see  (2.10)). 

The  integral  (2.3)  shows  two  saddle-points  at  ot  ■  -  it . 
For  an  asymptotic  evaluation  of  the  integral  for  kp  — »a> 
is  therefore  convenient  to  deform  the  original  integration 
contour  r  into  the  two-branch  steepest-descent  path  y 
(see  Fig. 2)  : 


=  i 


-  sg(a  ")  cos-1  ( 


cosh  a  " 


) 


(2.9) 


till'll 


/:  i ; 


i  in(a+<{ )  -si  nt}> 


iHu+i}> ) 


exp(-iki  cosalda 


♦  2  (Residues  of  F(a+4>)  in  I'-y  ) 


( 2 . :oj 


feq.  (2.10)  provides  the  ( Muliuzhinet s)  solution  to  the  pro¬ 
blem  of  the  plane  wave  normal  incidence  of  the  scattering 
geometry  depicted  in  l-’ig.l.  Discussion  about  the  location 
of  poles,  and  evaluation  of  corresponding  residues,  is  de¬ 
ferred  to  Appendix  B.  Hereafter  only  conclusions  are  sum¬ 
marized,  with  reference  to  radiowave  propagation  over  the 
l.arth,  so  that  II  polarized  waves  are  of  interest,  and 

1 

s  i  nO,  =  - -  ,  _ 

h  /;  +;w,V”"  (2.11) 

» G  +  10  /  Wt 

*  0 

wherein  e  and  j  are  the  appropriate  relative  permitti¬ 
vity  and  conductivity  (different  for  y  >  0  and  y  <  0) 

of  the  liarth. 

With  reference  to  Fig. 3  (wherein  the  angle  4>q  is  assumed 
real,  for  simplicity  J,  residues  of  F(a+<|>  )  correspond 

to  the  incident  wave 


Fi(a»0)  *  exp(-ikpcos(4>-4»0)) 

the  reflected  wave  for  y  >0 

+ 

r  cos 0  -sine 

V  (h  ,  p  )  =  - - _ ,.i 

+  +  F  exp(jkpcos(<f+  ♦,  )  ) 

costf  +s  inu  o 

existing  only  for  r  >  the  reflected  wave  lor  y  >0 


Fr(  r-  ,d>)  = 


:os$0  -  sin  0 


cos<{>  +sin  3 


1'0  expUkpcosUH  )) 


(2.12) 


(2.13) 


(2.14) 


existing  only  for  f  .The  integrand,  however,  posses¬ 

ses  other  poles,  which  reside  outside  the  closed  contour 
r-y  but  can  approach  the  saddle  points  at  1  it  so  to  in¬ 
fluence  the  asymptotic  evaluation  of  the  integral  in  (2.10). 
This  evaluation  is  deferred  to  in  Sect.  4,  after  the  solu¬ 
tion  (2.10)  is  extended  to  the  oblique  incidence  case. 
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3.  Scattering  under  oblique  incidence.  Extension  of 
Ma liuzhinets  'theory . 

Let  us  now  extend  the  two-dimensional  solution  discussed 
under  Sect. 2  to  the  three  dimensional  problem  of  oblique 
incidence,  wherein  the  incident  wave  is  a  plane  wave  of 
components  of  the  type: 

F(p,<(>,z)  *  F*  exp(ikzcosB)  exp(-ikpcos(<M>0))  (3.1) 

depicted  in  Fig. 4.  Without  loss  of  generality  we  will  as¬ 
sume:  0<  Re(B)  _<  rr / 2  ,  0<  Re(4>Q)  <  rr/2  and  will  drop  the 
z-dependence  for  the  field  components. 

Boundary  conditions  are  now  expressed  as 

L-xx*  i!  *  Z±xxn  at  <j>  =  ±u/2  (3.2) 

or,  equivalently, 

H-X5MI  «  (l/Z±)Lxx  at  ♦-tir/2  (3.3) 


By  taking  the  divergence  of  (3.2  and  3)  and  using  Maxwell's 
equations,  we  get  : 


-  0  at  <j>»±7r/2 

P  3<J>  h 

1  3Hv 

- -riksin0^Hx  *  0  at  4>"±7r/2 

p  9(j) 


(3.4) 


respectively,  fiq.s  (3.4  and  5)  are  symilar  to  (2.1)  and  sug¬ 
gest  that  results  of  Sect.  2  could  possibly  be  extended  to 
the  oblique  incidence  case  by  operating  upon  the  x-components 
of  the  field.  It  will  be  shown  in  the  following  that  this 
extension  is  not  straightforward. 

Let : 

' 

E  (o,4>)  =  — —  \  E  (a+4>) exp ( - ikosingcosn) da 
x  2-ni  /  x 

«  r  (3.5) 

Hx(o,$)  -  — j  H^a+tfOexpf-ikosinecosoOda 
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Then,  forcing  boundary  conditions  (3.4)  we  get  as  in  Sect. 2 


+  ~ 


(sina  -sin  e"  lE^fa*  tt/2)  *  (-sina±sin0h  )Ex(-a-u/2) 

i  -  ±  ~ 

(sina  -sin  6^  )Hx(a-  11 /2)  =  (-sina-sinO^  )  |j((-a-'n/2 ) 


(  3.6a  ) 
(  3.6b  ) 


wherein,  at  variance  of  (2.2), 

± 

+  Z  + 


sin  6, 


£sin  6 


sin  0 


sinB 


(  5.7  ) 


The  general  solution  of  (3.6a  or  b)under  the  condition 
that  a  single  pole  exists  at  a  *  <t>  in  the  s^rip  -tt/ 2  < 

<  a  <  tt/2  is  the  same  as  (2.6).  However,  and  11^  are 
no  more  necessarily  regular  for  p  - 0  ,  so  we  cannot  ap¬ 

ply  the  arguments  of  Sect. 2  .  It  is  therefore  convenient 
to  express  E,  , 1 1 2  in  terms  of  E  ,  llx  (for  details, 
see  Appendix  “  C)  ,  hence: 


f  ~  cosBcosaE  fa)  -  sina  £11  (a) 

E,(a)  ■  - 


l  dlz(a)  =  • 


1  -  cos2asinzB 

A  A 

cosBcosaHx (a)  ♦  sina  Ex(a) 


1  -  coszasinzB 


sing 


sing 


(3-8) 


Now,  (E  »H  )  should  approach  u  constant  value,  as|a"| 

— >  CD  . 


Ex(a) 


£Hx(a)  - 


Accordingly,  from  (2.6): 

r*  i 

Hox  cos  <t>n  I  i|»h  (°0 


wherein  E 


sina-sin$Q 
Wox  COE*o 


sina-sin<t> 


H 


♦  Bi. 


W 


<l'e(a) 


J  \j;  («)>  ) 


13.9) 


,  B,  and  B  are  constants,  and  (a)* 
ox  ox  he  .  zn 

^(a.0h)  ,  <Pe  =  <Ku,6e)  (see  Appendix  A).  (Ejx,Ei  )  ard 

the  values  of  the  x-components of  the  incident  field  at  p»  0, 

and  a»*e  easily  computed  forcing  the  residues  of  the  spectral 

representation  for  (E2,llz)  at p  “0  to  coincide  with  the 

values  of  the  incident  components  (Eot,Hoi^  at  p  hence: 


k 


i 
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'  l.1 

ox 


^02  cosgcos»0  +  CH*Z  sin4>0 
sing 


w 

V  ox 


r)\Qz  cosBcos<J>o 


E  sin4> 

02  O 


sing 


(3.10) 


For  evaluating  the  remaining  constants  (B^B^)  we  note  that 
the  denominator  of  (3.8)  vanishes  for 


2 

cos  a  = 


(sin2g)_1 


(3.11) 


For  g  real  eq. (3.11)  is  verified  for 


nn 


-  i  cosh  *  ( 


sing 


-) 


(3.12) 


Solutions  of(3.11)for  complex  g  are  provided  in  Appendix  D, 
but  they  are  not  relevant  to  the  discussion  which  follows. 


When  the  integration  path  T  of  the  integrals  (3.5)  is  de- 
formed^onto  the  steepest  descent  path  y,  the  two  poles  at 
a  a  a1  -  <j>  are  crossed.  The  corresponding  residues  will  describe 
(inhomogeneous)  plane  waves  of  the  type 


exp(-ikpcos4>)  exp(  tkpsin^cosg) 


(3.13) 


It  is  evident  that  such  two  waves  are  diverging  in  f  i  o 

o 

and  therefore  should  be  cancelled  if  the  solution  is  required 
to  be  bounded  everywhere  in  space.  This  can  be  accomplished  by 
forcing  the  numerators  of  (3.8)  to  vanish  for  a  “a*  ,  hence: 


Ev(a*  )  *  itH  fa*  )  -  0 


(3.14) 


In  deriving  (3.14)  note  that 

sin  a*  -  *  i  cotg  8  (3.15) 

as  it  follows  (3.11) 

Eqs  (3.14)  are  a  system  of  two  equations  in  the  two  unknowns 
(  B  8.)  and  the  solution  is: 

Of  h 
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B  = 
e 


eosi{>0  sing 


1  -  sin'gcos^  <J> 
cos4>0  sing 
9  2 

1  -  sin^gcos  <p 


(sin4>  sing-  iQcosg)-  II ‘.S  cosg 

vA  U  Ua  vT 


(3.19) 


(sin4>Qsing+  iQcosg)  *  E^S^osg 


wherein 


Q=  - 


S  = 
e 


^e(ao^h(ao^  +  ♦e^o^^o5 

2  <l>e(a*He(a') _ 

+  ^e(aA(ao) 

„  = _ 2  _ 

il  (a  )tfv(a  )  +  (a 

ye  o  yh  o  ye  o  h  oJ 

4.  Asymptotic  evaluation  of  the  field. 


(3.17) 


(3.18) 


(3.19) 


With  reference  to  the  two-dimensional  case  considered  in 
Sect.  2,  the  explicit  evaluation  of  the  field  requires  com¬ 
putation  of  the  integral  which  appears  in  (2.i0).  When  kyo  is 
large,  asymptotic  techniques  can  be  applied. 

In  order  to  have  a  uniform  asymptotic  evaluation  of  the 
integral,  ca*a  should,  be  taken  concerning  the  possible  loca¬ 
tion  of  poles  of  the  integrand  nearby  the  two  saddle  points 

a  *■  f  it 
s 

Poles  of  the  integrand  do  appear  at  (see  Appendix  B)  : 

<*n  =  -  <t>  +  (~l)n4>0  +  nir  ,  n  *  0;-l;i2; - 

± 

a+  =  -  4>  ±  (3tt/2  +  0  ) 

It  is  evident  that  and  _1T  when  $■*•<!> 0 

Furthermore,  when  the  incidence  is  at  grazing  angle,  i.e., 

<j>  'vti/2.  dnd  the  observation  point  is  close  to  the  surface 
y<0.  ija.,  4>  ^  v/2.,  then  a0-*-  n  .  • 

Under  these  circunstances ,  and  e~+  0,  also  a_  +-v  . 
Accordingly,  the  conclusion  i*  drawn  that,  #©*  H  -  polarized 
waves  and  radiowave  propagation  over  the  Earth  surface,  the 
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poles  «o,aj  approach  the  saddle  point  as  =  tt  ,  while  the 
poles  ot_2  ,  a_  cluster  around  the  saiuie  point  ots  =tt  . 

On  the  other  hand,  this  is  a  very  important  situation 
for  practical  application,  when  both  the  source  and  the 
observation  points  are  close  to  the  iiarth  surface,  for  y>0 
and  y< 0  respectively . 

Accordingly,  we  will  present  in  this  section  an  asymptotic 
evaluation  of  the  integral  which  appears  in  ( 2 . 1 0)  for  an  H- 
polarized  incident  field  and  tt/ 2  ,  0'v -tt/ 2  ,  0~<<  1. 

Let  us  first  consider  the  integral 


j  - 

J  s  in( 


os  4> 


(u+:J))-sin4> 


iKot+4>) 


<H4>q) 


exp(-  ikt>eosa)da 


(4.1) 


whose  integration  contour  ^  ,  is  the  branch  of  y  passing 
through  the  saddle  point  a  =  a,  =  it  ,  with  nearby  singularities 
at  a  =  aQ  and  a  =  a  j  . 

A  simple  procedure  for  the  asymptotic  evaluation  of  the  in¬ 
tegral  (4.1)  will  now  be  presented;  a  more  complete  one,  a  mo¬ 
dified  version  of  Bleistein  method  jsl]  (  which  allows  the 
error  estimate) is  given  in  Appendix  F. 

Let  us  expand 

cos*  -J  sina  ]  sina, 

- 2. -  .  _J - S__,  ■■  —  -  .  B«.)  0-2) 

sin(a+<)>)-sin<}>o  cos4'ia-cos4']a0  cos“lu-cos“iaj 

Wherein  B  (  a  )  is  regular  in  the  neighbour  hood  of  a=iT  , 

When  (4.2)  is  substituted  in  (4.1)  the  main  cent ribution  to 
the  integral  is  provided  by  the  first  factor,  when  <f>  n/2  , 

4>o  %  tt/2  ;  therefore,  we  can  neglect  the  contribution  arising 
from  the  B(u  )  term.  Then  we  are  reduced  to  consider  the  ca¬ 
nonical  integral 

]  sin  a  (|ra+<J>) 

- - - - - exp(- ikpcosakla  r  j  -j . 

:os“ ja-eos  jo  C«t>0 ) 


whose  solution  is  given  in  Appendix  H. 

The  total  field  H(  P  ,  $)  is  then  represented  as  a  sum  of  five 
terms,  which  will  be  discussed  hereafter. 

The  first  term  is  the  sum  of  the  direct  field  (2.12)  and 
the  contribution  of  the  integral  (4.11)  at  a  *  ->)>+4'0  : 

;  i  F(-/zto  cos]  (*  -ifd  ) 

H  U’.d)  -  H  exp (-ikpeos (<{>  -()))) - 

0  0  2F(0) 


■  II*  exp ( ikpeos (fio+6)) 


F(-/2kp  sen]  ( tfo+6) ) 
2F(0) 


(4.4) 
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where  the  two  convenient  (.positive)  angles  *  /■ 


*  =  n/2-6  ;  4>  =-n/2+6 

o  o 


(4.5) 


have  been  introduced, The  field  constituent  (4.4)  is  the  in¬ 
cident  wave  (2.42)  "corrected"  by  the  (normalized)  FrSsnel 
integral,  and  is  symbolically  sketched  in  the  first  line  of 
fig.  5. 

The  second  term  is  the  sum  of  the  reflected  field  (2.43) 
at  y  > 0  and  the  contribution  of  the  integral  (4.1)  at  a  ^ 


=-(4>0  +  <J>  ) 

H*(i  ,  ..) 


+if : 

+ 

sin«S  -sin  0 
o _ 

sin5  +sin  6 
o 


exp(ikpcos(60- 6)) 


F(/IEp  sinU^Q-^) 
2  F  ( 0 ) 


(4.6) 


Where  6+  is  a  shorthand  notation  for  e  . 

The  field  constituent  (4.6)  is  the  reflected  wave  (2.13)  with 
the  appropriate  (normalized)  Fresnel  transition  function  which 
renders  this  field  continents  for -V 2<4><tt/ 2  and  is  symbolically 
sketched  in  the  second  line  of  Fig.  5. 


The  third  term  is  the  sum  of  the  reflector  field  (2.44)  at 
y  <  0  and  the  contribution  of  the  integral  (4.x)  at  a  = 

—  7T  —  f  fK  4-  A  “1 


‘-n-(4>  +  <fo) 
H"(p,d> )  = 


sino  -sin  e" 


sin6o+sin  0 


-  H  exp(ikpcos(5  -6) ) 


F(/2l«r  sin](<S-6  )  (4.7) 


2F(0) 


The  field  constituent  (4.7)  is  the  reflectod  wave  (2.14)  with 
the  apnropriate  (normalised)  FrSsnel  transition  function 
which  renders  this  field  Continous  for  -tt/2<4><  it  /2  , 

and  is  symbolically  sketched  in  the  third  line~of  Fig.  5. 

The  fourth  term  is  the  contribution  of  the  integral  (4.1) 
at  a  -  -  4* ■* 3 7f / 2  -eT  : 


H_(p,<p) 


2  sin<5  tg  0  ^(tt/2+0  ) 


cos6  -COS0 


iP(tt/2-6q) 


Hq  exp(ikpcos(6+0^) 


F(v-^te  sinj (<5+0^) 


2F(0) 


(4.8) 


The  field  constituent  (4.8)  is  the  surface  wave  excited  on  the 
half-plane  y<  o  "  corrected"  by  the  (  normalized)  FrSsnel  in¬ 
tegral  and  symbolically  sketched  in  the  fourth  line  of  Fig.  5. 
The  fifth  term  is  the  contribution  of  the  integral  (4.1)  at 

a_2  =  +4>q  ~  2t r  ' 


A 
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He(p,<j>) 


sin6_  -  sin6 


sin6  +sin  6 
o 


sin6o+  sine 


sin6  -  sin  6 
o 


—  H*  exp(ikp  cos(6-*-6q)) 


F(^Fp  sinl(60+6i)) 
2F(0) 


(4.9) 


The  field  constituent  (4. 9) shows  the  same  symbolic  repre¬ 
sentation  as  the  field  constituent  (4.4),  i.e.,  the  first  line 
of  Fig.  5.  In  order  to  understand  its  physical  leaning,  Jet 
us  sum  and  substract  the  incident  field  (2.12)  to  the  total 
field.  Then. 

i  i  sin} (6„-6)) 

H(p3  «*>)  *  tr  exp ( ikpcos ( 6q+6 ) )  +  R+H*  exp(ikpcos(6Q-6))  - 


+  R_H^  exp( ikpcos (6-6  )) 


2  sin6  tg6 
o 

cos6Q-cose 


o 

iKit/2+6  1 
iP(tt/2-60) 


F(/5kp  sin}(6-60)) 
2F(0) 


H*  exp(ikp  cos(6+e  )) 


2F(0) 


F(/2kp  sin} (6+0)) 


2  sin60(sine'*-sine+) 


(sin6Q*sine+) (sin60-8ine“) 


Hj  exp(ikp  Cos (6  -6)) 


2F(0) 

F(v^kp  Sln}(6^6))  (4  10) 

2F(0) 


Accordingly,  the  total  field  for  4»^-tr/2  ,  ♦0^it/2  is  the  sum 
of  the  incident,  the  two  reflected, (  with  appropriate  transi¬ 
tion  function)  the  surface  field  and  an  extra  (pole-type)  field 
due  to  the  discontinuity  of  the  surface  impedance  at  y=0 
In  general, also  the  diffracted  wave  (saddle-point  type)  due  to 
the  discontinuity  should  be  considered  (see  Appendix  F) . 

5 .  Field  produced  by  a  line  source  over  a  mixed-path  plane  Earth. 


Let  us  nOW  consider  an  uniform  magnetic  line  source  parallel 
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to  the  surface  impedance  discontinuity,  as  depicted  in  Fig.  6  , 
and  located  at  the  point  ty[x  ,  y  ). 

The  tree-space  magnetic  field  radiated  by  the  source  in 
F ( \ , y )  would  be  given  by 


flV.y)  -  H1(p,«>)  - - - —  Hll\kd)  -  (kd) 

4c  o  0 

wherein  lm  is  the  magnetic  current. 

By  using  the  integral  representation  for  the 
we  have 


HV,<S) 


exp(ikd  cos (6  -n))d6 
o  o 


llanke  1 


(5.1) 


func t ion 


—  /  exp(ikp  cos(6  -n  ))exp(ikp  cos(6+6) )d6  (5.2) 

■p  /  u  u  u  O  0 

^  S 


wherein  r  is  the  Sommer feld  integration  contour  as  depicted 

in  Fig.  7  and  we  have  assumed,  for  convenience,  y  >  y. 

‘  o  ' 


Eq.  (5.2)  represents  the  field  radiated  by  the  line  source 
as  a  superposition  of  plane  waves  with  spectral  distribution 
exp(ikpocos(60  -n0)) -Accordingly, the  field  H*  radiated  by 
the  source  in  the  presence  of  the  mixed  path  configuration 
of  hig.6  is  given, by  superposition  , as 


Hs(p,6) 


C  f 

m  -  /  H(n, 5,  5Q)exp(ikp0  cos(6o-n0))d6(. 


(5.3) 


wherein  H(p»5»5p)is  given  by  (4.10''  plus  the  possible  con¬ 
tribution  from  the  field  scattered  by  the  discontinuity  at 

y-o. 

I  he  perturbation  in  the  field  excitation  due  to  the  impe- 


a 
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dance  discontinuity  is  of  interest. The  change  in  the  spec¬ 
trum  can  easily  be  obtained  by  assuming  6  =9  (the  surface 
impedance  is  now  costant  and  equal  to  V  ever ywher ej in  (4.10); 
and  then  subtracting  resulting  expression  from  (4 . 10)  , hence : 

sin9  -sin6+  .  fexp(ikpcos(6  -6))  F(/2kp  sin£(6  -6)) 

All  =  2  sin6  - --  ^ - 5 - 

sin<5o+sin9  |sin6o+sin9  2F(0) 


exp(ikpcos(6Q+6))  F(»^kpsinJ(<s0+|5) ) 


(5.4) 


sin6  -sine 
o 


2F(0) 


For  6  =0  (observation  point  upon  the  Earth's  surface), 
expression  (5.4)  further  simplifies  as 


4sin  6Q(sin0  -sin6  )  .  F(/2kpsinK60' <$)) 

ah  - ? - 2 — rr^r  Ho  exp(ikpcos6o}  — 

(sin6Q+sin6  ) (sin  <$o-sm  8  )  2F(0)  ^  ^ 

Since  p  ccs(6  -n  )  +  pcos60  =dcos  (60  -n)  ,we  have  the  integral 
representation  for  AHS  . 


4C(sin0~-sinO+) 


F(/2kpsinl6Q) 

2F(0) 


sin26_ 


(sin60+sine+) (sin26o-sin2e”) 


exp(ikd  cos(6Q-n))d6o 


(5.6) 


i.e.  for  the  perturbation  of  the  magnetic  field  on  the  Earth 
surface  due  to  a  discontinuity  in  the  surface  impedance  and 
excited  by  a  magnetic  line  source. 

Analytical  evaluation  of  the  integral  (5.£)|s  not  an  ea¬ 
sy  task. The  integrand  presents  poles  at  6  =-8  ;-9  , clustering 
around  the  saddle  point  <5  =n  (when  n<<l)TSis  complication  can 
be  overcome  using  the  evaJuation  techniques  presented  in  Sect. 

4 .However , the  integrand  is  also  rapidly  varying  due  to  the  pre¬ 
sence  of  the  Fresnel  integral  ,when  kp»>l.This  point  requires 
further  study  to  be  overcome. 
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6 .  Conclusions  and  recommendations 

The  problem  of  plane-wave  scattering  by  a  surface  impe¬ 
dance  discontinuity  over  a  plane  has  been  rigorously  solved. 

The  solution  is  attractive  inasmuch  it  provides  a  physical 
description  of  the  plane-wave  scattered  constituents.  These 
consist  of  the  pole-type  contributions,  essentially  plane  wa¬ 
ves  with  attached  Fresnel  transition  functions;  and  saddle- 

point  type  contribution,  essentially  a  cylindrical  wave  scat¬ 
tered  by  the  discontinuity.  For  an  observation  point  close  to 
the  surface  and  an  incident  wave  at  a  grazing  angle,  a  novel 
result  is  the  appearance  of  a  new  type  of  (  pole-type  )  con¬ 
tribution,  in  addition  to  the  conventional  incident,  reflected 
and  surface  fields.  The  role  of  the  (  cylindrical  )  wave  scat¬ 
tered  by  the  discontinuity  should  be  better  analysed,  and  graphs 
pertinents  to  a  pratical  situation  computed. 

The  field  excited  by  prescribed  sources  can  be  computed, 
in  principle,  from  the  plane-wave  expansion  of  the  incident  field 
and  then  application  of  results  of  the  scattering  analysis.  How¬ 
ever,  the  resulting  integral  representation  contains  a  rapidly 
varyng  term  (  Fresnel  integral  )  in  the  neighbourhood  of  the  sad¬ 
dle  point,  and  also  as  many  as  three  poles.  This  is  a  point  which 
needs  further  attention  and  the  use  of  appropriate  (  if  available  ) 
mathematical  techniques  of  asymptotic  evaluation  of  integrals. 


Appendix  A  :  The  function  ip  Coi)  . 


The  function  tp  (ex)  is  defined  us  follows  [20  ) 


4>a)  -  N  ( a  n  -  0  )  N  (  a  ♦  0  1  N  (.  a  -  tt  +  6  )  N  (  a  -  0  ~  ) 


1 A  .  1 


*r  equivalently: 


t  (oil  =  N  “  |  u  /  ’ 


N(a  +  0  )N(a-6  1  a_0  +  +  11/2 

—  c  o  s  - 


ot  +  P  -  it  / 


c  o  s 


N ( a • 9  )N(a+0* ) 
The  function  N  (.  ot  )  is  given  by 


(,  A .  2  ) 


a 

1  r  2t-Trsin  t 

N(a)  =  exp  (  -  - dt  )  (/.3) 

4  TT  J  COS  t 

0 

and  is  a  particular  determination  of  a  most  general  class 
of  function  introduced  by  Maliiuhinets  (20]  with  reference 
to  the  problem  of  impedance  wedge  scattering.  The  function 
N  (  a)  is  easy  to  compute  numerically  and,  in  any  case,  it 
is  simply  related  to  a  tabulated  function  [21,28].  Its  main 
properties  are  hereafter  summarized: 


N  ( a )  =  N  ( -  a )  ( A .  4  ) 

N  (  a  tt  /  2  )  Nla-n/21  =  N 1  (  tt  /  2  )  cosa/ 2  S) 

1  im  , 

ju"|  *ocN(a)  ~  0  (  exp|a"|/4)  fA.b) 

The  location  of  poles  of  i{((a)  within  the  strip 
-  3  tt  /  2  <a'<  3  tt  /  2 


is  of  interest.  From  (A. 3)  it  is  notod  that 
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2t  -  usin  t 

-  dt  (A. 7) 

cos  t 

is  finite  for  a  =  ir/ 2  ,  diverges  toward  negative  values 
for  a  =  3it/2  ana  toward  positive  values  for  a  *  5 rr/ 2 
Since  the  divergence  of  (A. 7)  is  of  logarithmic  type, 
the  function  N(a)  has  a  simple  pole  at  a=±  Stt/2.  Accord 
ingly,  the  poles  of  ip(a)  in  the  strip  of  interest  are, 
from  (A.i): 

a±  =±(3ir/2  +  6‘)  (a.r) 


S 


0 


Appendix  Bj_  Poles  of  the  spectrum  of  the  field  and  re¬ 
sidues  of  the  integral  representation. 


Poles  of  the  trigonometric  factor  in  (2.8)  do  oc- 
cour  at 


otn=’  ~<t>  +(-l)n0o  +nir 


ns  0  ;  -1  ;  . . . .(B-l) 


For  n=0,  the  residue 

P1  (p,<|>)  =  F*  exp(-ikpcos($-<J>0))U(ir-$ 

describes  the  incident  field,  wherein: 


♦e0)U(ir*4>-i-0) 

(B.2) 


5o  ■  *(5  *  .  ♦  „  "  •  ♦"  >_0  (B.J) 

U(t)  is  the  Heaviside  unit  step  function  and  gd  (x)  is 
the  Gudermanian  function 

gd(x)  =  sn(x)  cos_1( - - - )  =  s in-1 ( tgh  x) 

cosh  x 

For  nE  -  1,  the  residues  (B.4) 

) 

fI  (p*4>)  ■  -  F*  - —  expCikpcos  ( «J>*4>  ) )  • 

1  °  ° 

■U(2ttT(C0  +  4>))U(*(C0*(D)) 


CB  .5) 


describe  the  reflected  fields  upon  the  two  half-planes 
t  n/ J  respectively.  Use  of  (J.5i  shows  that: 


o,  ' i tt - 4>  )  cos4>  -sine  i 

-  — - - —  = - ° - - —  *  R  (n.ft) 

^  *  ■  j  *  cos^*  s  i  nO 

Accordingly,  the  excitation  factor  for  the  reflected 
waves  is  identified  with  the  Fresnel  reflection  coeffi¬ 
cient  R*  appropriate  to  the  two  sur  laces  <p  v/2>  as 

it  should  be  expected.  For  an  incident  ll-polar  i  zed  sur¬ 
face  wave  at  the  Brewster  angle,  i.e.,  4> o=  v/2- 
there  is  no  reflected  wave  at  #  =  tt  /-  ;  the  reflection 

coefficient  at  <J>  =  _rT/2  takes  the  simple  form 


1‘olcs  of  the 


+ 


sin0h  -  sin0h 

+ 

sin0h  +  sin6h 
other  factor  in  (2 
a*  *  -  $  *(371/2  + 


.8) 

+ 

0" 


do 

) 


IB.?) 


occour  at 

( B  .  8 ) 


as  follows  from  (A. 7).  The  corresponding  residues 


F  (  p  . 


2cos4iotg0 

sind>  *cos  0 
o 


^[+(0~  +  Tl/ZlJ 


exp  [*  ikp  s i n C T  8  )J  U ( -  C"  £  ~  "  r  /  2 ) 

(  a  •  d ) 

5  =  e '  +  gd(e”3  ,  0=  o '  +  i  o"  ,  e  ’>  o 

correspond  to  surface  waves,  one  (upper  sign)  with  phase 
velocity  in  the  direction  tt/2  +  0+  ;  and  the  other  (lower 
sign)  with  phase  velocity  in  the  direction -n/~  -0“ 

In  computing  (B.9)  use  has  been  made  of  (2.5).  These 
surface  wuves  do  exist  in  the  angular  range 

>  7T/2+  5* 


(  B .  1 1  / 


Ihen,  a  necessary  condition  for  such  waves  to  exist  is  that 


S  =  9  '  +  gd(0")  <  0 

Letting  sine  =  r  *  i;<  ,  tedious  algebra  shows  that  this  is 
is  equivalent  to  the  condition  : 


-  x 


f  B  .  i  2 ) 


For  radiowave  propagation  over  the  liarth,  li  -  polarized 

waves  are  of  interest  and  1 

Sin  9h  =  7  ~T~C  (rt.13) 

03CO 

wherein  c  undo  are  the  relative  permitivity  and  conducti¬ 
vity  of  the  Larth.  Then  x  <  o  and  (B.12J  transforms  as  fol¬ 
lows  : 

2  -i 

)  (B .  14) 

o 


m 


The  conclusion  is  drawn  that  surface  waves  (B.9)  are  not  excited 
since  the  poles  (B.8)  are  not  crossed  upon  deformation  of  the 
intergration  contour  from  r  to  y 


Appendix  C:  Spectral  relations  between  longitudinal  and 


transverse  field  components. 

'  _  expfikz  cosB)  r 

:  z  =  —  J  Ez(a+0)  exp(-ikpsinBcosa)da 

j  "  r 

1  (C.l) 

exp ( ikz  cosg)  r 

z  p,’*>,z-1  =  ~  J  flz(a+tf)  expf-ikpsingcosalda 

r 

It  is  known  (30)  that  the  total  field,  (Ij  ,  H  ),  can  be 
expressed  in  terms  of  (  L  ,  H_)  : 
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..  .  -2  .  -2 

h  =  k  sin  g  (VxVxE  z  +  icon  xH  z  ) 

z  z  ' 

(0.2) 

-  2  -  2 

H  =  k  “sin  g  (vxvxfi  z-  i ui c  xE  z  ) 

z  z  ' 


Substituting  (C,l)  in  (3.2)  ,  projecting  on  the 

x  -  axis  and  equating  the  corresponding  spectral  com¬ 
ponents  ,  we  get 


Ex(a)  = 

?Hx(a)  = 

The  system 
hence  : 


JS  A 

cosg  cosa  E  (a)  +  sina  (a) 
sing 

A  A 

cosgcosa  ^Hz(a)  -  sina  E„(a) 
sing 

(C.3)  can  now  be  solved  for 


(0.3) 


< 


E 

•7 


H 


z 


cosa  cosg  E  -  c  sina  Hv 

X _ X 

2  2 

1  -  cos  a  sin  g 
cosgcosa  CH  +  sina  £ 

_ A X 

,  2  .  2  „ 

1  -  cos  a  sin  g 


sing 

(C.4) 

sing 


.  Appendix  D  -  Spurious  poles  in  the  case  of  oblique 
incidence 


Poles  of  (3.11)  are  given  by 

cos  a  ■  i  (sing  ) _1 

Letting 

(sin  g  )  1  =  a  +  i  b 


(D.l) 

(D.2) 


Eq .  (0.1)  becomes 


{ 


cos  a' 
sin  a' 


cosh  a"  = 
sinh  a"  = 


4  a 
+  b 


(D.3) 


T 
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Where  F  (u)  is  the  FrSsnel  integral 


F(u)  = 


J ' exp(it2)dt 


-  1  according  to  0^  +  ga(0g)  T  0 


(1-4) 


(h.5) 


and  the  property 


F(u)  +  F  (-U)  »  2  F  (o)  */tT  exp  in/4  (l-h) 

has  been  used.  Note  that 

-f  oo 

/exp(-fit2)  2/n 

^—•3  dt  *  exp(-  ii—  )—  F(/5z)  (L .  7) 

—  oo 

The  other  integrals  can  be  computed  similarly. 


Appendix  F  :  Modified  Bleistein's  method  for  the  asymptotic 


evaluation  of  the  integral. 


Let  us  consider  the  function  B(a  )  defined  in(4.2)  : 


B(cO 


cos  4>0 

sin(a+4>)-sinci 


J  sin  a0  i  sin  aj 

cos^ Ja-cos^ laQ  cos2 Ja-cos^Jaj 

(F.l) 


and  let  us  expand  it  as  follows 


B (a)  ■  a  +  cos  }a  Bj  (a) 


(F.2) 


wherein  B(a)  is  not  only  regular  around  a  -  ag  *  n ,  as  B  (o  ) , 
but  also  attains  a  zero  value  for  a  ■  a  ■  n . 

Accordingly  in  the  asymptotic  evaluation  of  the  integral 


I 


2 


ip  (4+a) 
4(40) 


B(a)  exp  (-ikpcosa)  da 


CF.3) 


the  dominant  contribution  *  will  be  provided  by  the  first 
factor  at  the  right  hand  side  of  (F.2). 

Letting  a  =  tt  ,  we  can  easily  compute  a;  hence: 


-cos  4> 


a  = 


sin  4  +  sin  <p 


\  sin  (  4o-40 

cos2i(<J>0-<t>) 


1  sin  (<pQ-<p) 

cos2H<J>0+4») 


-  3  cosfto _ 

sin  1  (4q+4)cos  \  Oo-<J0 


then , 


I 


4(0+4) 

4C4J 


exp(-ikp  cos  a)  da 


I 


4Ctt+4+b) 


a  |  — -  exp(ikp  cosg)  dg 

4(4„) 


x  a 


■f  oo 


/2 


exp[i(kp-ir/4)J  i  'KB+n+<p)  exp(-kpt2) 

S0  ...  f  x  N  5 


-oo  4(40)  cos  j 


dt 


+  00 


1  ®  4  f  4+rcl  r> 

~1=~  exP  (.i  (kp-tr/4)J 

2  *(♦„)  J 


exp(-kpt2)dt 


*  xa 


4(40) 


(F.5) 
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This  corresponds  to  a  cylindrical  wave,  scattered 
discontinuity  at  y-0,  which  decays  with  distance  as 


-  3  cos  d> 

_ _ _  o  _ 

*'Tp'  sin]  C4>+«J>o)cosi  ( 4> +4>0 ) 

Accordingly,  its  contribution  to  the  total 
general  be  neglected. 


(F.6) 


field  c 
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Fig. 2.  Integration  contour  in  the  complex  a  plane. 


Symbolic  sketch  of  several  constituents 
field  close  to  the  surface. 


